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Abstract. Corrected trapezoidal rules are proved for J a f(x) dx under the assumption 
that /" E L p ([a,b]) for some 1 < p < oo. Such quadrature rules involve the trapezoidal 
rule modified by the addition of a term k[f'(a) — f'(b)]. The coefficient k in the quadrature 
formula is found that minimizes the error estimates. It is shown that when /' is merely 
assumed to be continuous then the optimal rule is the trapezoidal rule itself. In this 
case error estimates are in terms of the Alexiewicz norm. This includes the case when 
/" is integrable in the Henstock-Kurzweil sense or as a distribution. All error estimates 
are shown to be sharp for the given assumptions on /". It is shown how to make these 
formulas exact for all cubic polynomials /. Composite formulas are computed for uniform 
partitions. 



1. Introduction 

This paper is concerned with numerical integration schemes for J f(x) dx where it is 
assumed /' is absolutely continuous such that /" G L p ([a,b}) for some 1 < p < oo, or 
that /" is integrable in the Henstock-Kurzweil sense, or that /' is continuous so that /" 
exists as a distribution and is integrable using a distributional integral. Integration by 
parts shows that 

(1.1) J f(x) dx= l - [-f(a)4f(a) + f(b)<P'(b) + f'(a)cj>(a) - f{b)<t>{b)\ + E{f), 

where E(f) = (1/2) £ f '{x)<j>{x) dx and is a monic quadratic polynomial. Observe 
that taking <f>(x) = (x — a)(x — b) gives the usual trapezoidal rule J f(x)dx = (b — 
a) [/(o) + /(&)] /2 + E T (f), where E T {f) = (1/2) f a f"{x){x - a){x - b) dx. The Hdlder 
inequality then gives \E T (f)\ < ||/"|| p ||0|| g /2, where q is the conjugate exponent of p. 
Hence, 

r wrub-ay/s, P = i 

(1.2) |^ T (/)|<{ [B(q + l,q+l)f' ! \\r\\ p (b-af+^/2, 1< p < oo 

( ||r||oo(&-a) 3 /12, p = oo. 
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Here, B(x,y) = T(x)T(y)/T(x + y) is the beta function. See [U Theorem 3.22]. (This 
corrects a typographical error in [3].) 

We find (p that minimizes the error in (11.11) . This leads to a quadrature rule that includes 
values of / and /' at the endpoints a and b. The error-minimizing polynomial produces 
the classical trapezoidal rule, modified by the addition of first derivative terms. In the 
literature this is a called a corrected trapezoidal rule. This includes solving the problem 
of choosing k G R to minimize all quadrature rules of the form {b — a) [/(a) + fib)] /2 + 
k[f'(a) — fib)]. The error terms are as in ( 11. 2ft but with the coefficient of ||/"|| p (& — 
a) 2+l l q minimized. In particular, the coefficients are strictly less than in ( 11. 2p . We prove 
our results are the best possible given the assumption /" 6 L P i[a,b]) (Theorem 12.11 and 
Corollaries 12.21123]) . The composite formula (Corollary 12. 5p provides an improved error 
estimate over the trapezoidal rule. Since the /' terms telescope, the correction terms only 
require computation of /' at endpoints a and b rather than at interior nodes. Compared 
to the usual trapezoidal rule the extra computing time to implement our corrected rule is 
then negligible for large n. 

Finding the polynomial <p that minimizes the error in (11.11) involves solving a tran- 
scendental equation for a parameter in the polynomial. This transcendental equation is 
written in various ways in Section [3J When q is an integer this reduces to a polynomial 
equation for the parameter. We are able to solve for the exact value of the parameter when 
p = l,2,4/3,oo. See Corollaries I2.21l2~41 and ( 13. 3p . Even without knowing the parameter 
exactly, Theorem 13.21 gives a corrected trapezoidal rule with error estimate smaller than 
in flT2D. 

In Section [5] we reduce the assumption on / to / e C 1 ([a,6]) and then /" exists as a 
distribution and is integrable using the continuous primitive integral ( Corollary 15. 3j) . This 
includes the case when /" is integrable in the Henstock-Kurzweil sense (Theorem 15.11) . 
The error estimate is then in terms of the Alexiewicz norm of /". In this case, the optimum 
form of (II. ip is the trapezoidal rule itself. 

In Section[H]we compute <fr so that (11.11) is exact for all cubic polynomials /. The required 
<p is the same as the one that minimizes the error in the case when /" e L 2 i[a, b]). 

Several authors have considered corrected trapezoidal rules under the assumption that 
the derivatives of / are in various function spaces. See [I] (£ p ), [8] {L°°), [9] (Lips- 
chitz, continuous and bounded variation, L p ), [10] (L p , Henstock-Kurzweil integrable), 
[H] (continuous and bounded variation) and [17] iL p ). 



2. Main theorem 

The error in (II. ip is minimized over all monic polynomials <fi. Results in Lemma H~T1 show 
that a unique error-minimizing polynomial exists and is of the form 4>{x) = (x — c) 2 — a 2 
such that c is the midpoint of [a, b] and <j) has two real roots in [a, b]. Note that (II .ip 
becomes a corrected trapezoidal rule precisely when <p'ib) = —<f/(a) = b — a and this 
relation holds for 0(x) = {x — c) 2 — a 2 . For a uniform partition the composite rule 
obtained from (11. ip will in general have /' evaluated at all points at which / is evaluated. 
However, when 4>'(a) = 4>'(b) the sum of /' terms telescopes, leaving only f'ia) and fib) 
(Corollary 12. 5p . This is the case with all the error-minimizing rules we present. 
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We are able to compute exact values for a when p = l,2,4/3,oo. In other cases a is 
given by the transcendental equation (12.101) . When q is an integer this reduces to finding 
the largest real root of a polynomial of degree 2q — 1. 

Theorem 2.1. Let c be the midpoint of [a, b]. Let f: [a, b] — > M such that f is absolutely 
continuous and f" G L p ([a,b]) for some 1 < p < oo. Let 1/p + 1/q = 1. Amongst all 
monic quadratic polynomials used to generate (II. ip . taking <p(x) = (x — c) 2 — a q gives the 
unique minimum for the error \E(f)\. The constant (3 q > 1 is the unique solution of the 

f\x 2 - If' 1 dx = l -B{q, 1/2) = 2^ 2 B(q, q) 
(b — a)/2. This gives the quadrature formula 



equation 




(2.1) 




and a q (3 q 


= (b-a)/2 


(2.2) 


/ f{x) dx = 

J a 


where 




(2.3) 





6-ci (b — a 



[/(«) + /(&)] + ^-g^(l - Pf)\f(°) - /'(&)] + E(J), 
lF(f)l \\f"Ub-a) 2+1 ^-Pf) 

I UJI - 2 3 + 1 /9(g+ 1/2)V</ 

T/ie coefficient of ||/"|| p £/ie error bound is the best possible. 

Note that the numbers a g and (3 q are independent of /, while f3 q are also independent 
of the interval [a, b] . 

Corollary 2.2. If p = 1 and q = oo then (5^ = a/2 , aioo = (6 — a)/{2y/T) gives the 
unique minimum error. The quadrature formula is 

(2-4) T /(*) = <t=A [/(a) + /(&)] + ^£ [f (a) - /'(&)] + J3(/), 



2 l« \ / « \ / j lg 
where the optimal error is \E(f)\ < — a) 2 /16. 

Corollary 2.3. If p = q = 2 then (3 2 — V^, a 2 = (b — a)/{2\fT) gives the unique 
minimum error. The quadrature formula is 

(2-5) f f{x) dx = [f{a) + f{b)\ + ^£ [/'(a) - /'(&)] + £(/), 



2 L " v ' ' WJ 12 
where the optimal error is \E(f)\ < \\f"\\2(b — a) 2 ' 5 /(12v / 5 _ '). 

Corollary 2.4. If p = oo and g = 1 i/ien /?i = 2, ol\ = (b — a)/4 gwes ine unique 
minimum error. The quadrature formula is 

(2.6) jf' /(*) ^ = [/(a) + /(&)] + W^!-\f(a) - f'(b)] + E(f), 

where the optimal error is \E(f)\ < ||/"||oo(& — a) 3 /32. 

Other authors have obtained corrected trapezoidal rules under the assumption /" G 
L p ([a,b}), generally with different coefficients of /'(a) — f'(b) in (12.21) and strictly larger 
coefficients of \\f"\\ P (b—a) 2+1 / g than in (12. 3p . See Cerone and Dragomir [4] equation (3.64), 
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where the coefficient of f'{a) — f'(b) is (b — a) 2 /8 for all values of p. In Theorem 3.24 for 
p = oo they have the coefficient of f'(a) — f'(b) as (b — a) 2 / 12. This coefficient is only 
sharp for p = 2 (Corollary El. The estimate \E(f)\ < [sup(/) -inf (f)](b-a) 3 /(24:y/E) is 
proved. Dedic, Matic and Pecaric [HI Corollaries 9, 12] also consider corrected trapezoidal 
rules with /" G L p ([a,b]). In their Corollary 9 they have the coefficient of f'(a) — fib) as 
(6 - a) 2 /12 and prove the larger estimate \E(f)\ < \\f"\\oo(P - a) 3 /(18y/3). The results 
of our Corollary 12.31 appear as their Corollary 12 (without sharpness). Similarly in [10.]. 

Corollary 2.5. For a uniform partition given by Xi = a + (b — a)i/n, < % < n, the 

composite formula is 

(2.7) 



f(x) dx = h —^- 
2n 

where 



n-1 



/(a) + 2 £/(*,) + /(&) 



+ ^T^ 1 ~ P?) \f{a) ~ f(b)) + E(f), 



\\f"Ub-a) 2 /(lQn 2 ), p=l 



i £(fl |<j Kv<oo 
||f |U(6-a) 3 /(32n 2 ), p = oo. 
The coefficient of ||/"|| p in the error bound is the best possible. 

Proof of Theorem \2.1[ Let be a monic quadratic polynomial. Integration by parts gives 
(II. ip . We are then led to minimize E(f) = (1/2) y ° f"(x)(j)(x) dx. By the Holder inequal- 
ity, \E(f)\ < (l/2)||/ // || p ||0||q. First consider a symmetric interval [—a, a]. By Lemma 14.11 
to minimize \\4>\\ q we need only consider (f>(x) = x 2 — a 2 for some a G [0, a]. Let = a/a. 
Define 
(2.8) 

G q (a) = \\<f)\\ q = 2 J \x 2 - a 2 \ q dx = 2a 2q+1 (1 - x 2 ) q dx + / ' (x 2 - l) q dx 
Note that 

(2.9) G' q (a) = Aqa 2q uf (1 - x 2 )^ 1 dx - / ' (x 2 - I) 9 " 1 dx 

We have G' q (0) = 0. Since the function a (-> j^ a (x 2 — 1) 9_1 dx decreases from positive 
infinity to zero as a increases from zero to a, it follows that G has a unique root in (0, a) 
and G q has a unique minimum at a q G (0, a). The first integral in (I2.9P can be evaluated 
in terms of beta and gamma functions. We have [16, 5.12.1] 



J (1 - x 2 )"- 1 dx = ~S(g, 1/2) = 2 2q - 2 B(q, 



, v t r(g) 



21X2 + 1/2) 

Let /3 g = a/a g . The required minimizing polynomial is then determined by the unique 
root (3 q G (1, oo) of the equation 

(2.10) f \x 2 -iy- l dx= l -B{q,l/2). 
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To compute \\4>\\ q , evaluate the final integral in (I2.8p . Integration by parts establishes 
the recurrence relation 

(2.11) / & - I)' * = ^1 -^hJ<S- I)-' <*■ 

Using this and the corresponding version with integrand (1 — x 2 ) q , we obtain G q (a q ) = 
2a(a 2 - a 2 q ) q /(2q +1). It then follows that 

I Fm > < ||ni P (2a) 2+ ^[l-K/a) 2 ] 

I 2 3 + 1 /?(g+ 1/2)V? 

Replacing a with (6 — a)/2 establishes the error estimate for interval [a, 6]. Using (f>(x) = 
(x — c) 2 — a 2 the formula (12. 2p now follows. 

With the Holder inequality, \$ a f"{x)4>{x) dx\ < \\f"\\ P U\\ q , the necessary and sufficient 
condition for equality is f"(x) = dsgn[0(x)]|0(x)| 1 ^ p_1 ^ ) for some d G R and almost all 
x G [a,b). See p31 p. 46]. Integrating gives f(x) = d £{x - t)sgn[(j){t)}\(j){t)\ l l^-^ dt, 
modulo a linear function. □ 

Proof of Corollary \2.2\ It suffices to consider the interval [—a, a]. When p = 1 write 
<p(x) = x 2 — a 2 ^. By Lemma T4. II we need only consider the case with two distinct roots in 
(—a, a), i.e., < a„o < a. We have 



Fco(4>) '■= II0IU = max|0(x)| = max(|0(O)|, 0(a)) = max(a oc ,a - a 

\x\<a 



It follows that ttoo = aj\p2. Formula ( 12. 4ft and the error estimate now follow easily. 
There is equality in |J" f"(x)<p(x) dx\ < ||/"||i||0||oo whenever <p(x) = dsga[f"(x)] for 
some del and almost all x G [—a, a]. See [131 P- 46]. This does not hold but we can 
show the coefficient of ||/"||i cannot be reduced by considering a sequence of functions. 
If /" = 5, the Dirac distribution supported at 0, then \f f"(x)(x 2 — a 2 ^) dx\ = a 2 ^ = 
a 2 / 2 = || /"|| i|| 0|| oo- Of course, if /" = 5 then /' = X(o,oo) which is not absolutely 
continuous. Let (i[) n ) be a delta sequence. This is a sequence of continuous functions 
4> n > with support in (0, 1/n) such that f^ n ' if) n (x) dx = 1. Now define f n (x) = 
J- J-a VV).^) dz dy. Then f n is absolutely continuous and ||/"||i = 1. Since is continuous 
we have lim n ^.oo| f"(x)(x 2 — a 2 ^) dx\ = a 2 /2 = ||0||oo- The error estimate is then 
optimal. □ 

Proof of Corollary \2.4\ Write 4>(x) = x 2 — a\. By Lemma [4.11 we need only consider the 
case with two distinct roots in (—a, a), i.e., < a± < a. Equation (12.91) now becomes 
G'i(ai) = 2ai(12ai — 6a)/3, from which a± = a/2. This then gives ( I2.6p . There is 
equality in \ f"(x)4>(x) dx\ < ||/"||oo||0||i whenever f"(x) = c?sgn[0(a;)] for some d G R 
and almost all x G [—a, a]. See [131 P- 46]. Integrating gives f(x) = d f* a (x— £)sgn[0(t)] dt, 
modulo a linear function. □ 

This case also appears in [221 Theorem 1]. 

The proof of Corollary 12.51 follows using the Holder inequality for series as in the proof 
of Theorem 3.26 in [4]. 

Lemma [4.11 shows that the minimum of ||0|| 9 over monic quadratics is unique. Hence, 
the coefficient of ||/"|| p (& — a) 2+l / q in (12. 3p is strictly less than for any other choice of 
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a q and indeed for any other choice of monic quadratic. In particular, we get a smaller 
coefficient than in the trapezoidal rule (II. 2p . 



3. Evaluation and approximation of /3, 



q 



In Corollaries 12 .21 through |2~41 we were able to compute the exact value of a q and /3 q for 
q = 1, 2, oo. In this section we compute /3 4 as the root of a cubic polynomial. When q is 
an integer, equation (12.101) becomes a polynomial of degree 2q — 1. See (13. ip and (13. 4p . 
When q is even the degree reduces to q — 1 and we compute the exact value of /3 4 (13. 3p . 
In general, equation (I2.10p is transcendental and most likely cannot be solved exactly. 
We rewrite this in terms of hypergeometric and associated Legendre functions and also 
show that v^T < (3 q < 2 and is decreasing (Proposition 13. ip . In Theorem 13.21 we show 
how the corrected trapezoidal rule can give good approximations of the integral of / with 
/" G L p ([a, b}) for all 1 < p < oo even if the exact value of f3 q is not known. Part (c) of 
this theorem shows that the corrected trapezoidal rule with a = 0, (3 = oo gives a smaller 
error estimate than (II. 2p for all 1 < q < oo. The coefficient is a simple function of q. 

If q is an integer, use the binomial theorem to write (I2.10p as 



fc=0 v 7 fc=0 v 7 

This gives a polynomial of degree 2q — 1 for /3 q . When q is even this reduces to 



1 7 ~ _ a (-i) fc ffi 

x fc J 2A; + 1 

fc=0 x 7 

?4 i 7fl2 



When p = 4/3 and q = 4 we get the polynomial /3| — (21/5)/3| + 7/3| — 7 = 0. The 
unique solution is [161 1.11 (in)] 

(2(7 1/3 ) fx , iV3 r , , iV3\ 1 1/2 

(3.3) 4 = | (32/3)5 ([5v / 3(r + 27j -pv / 3lT-27j J + 7/5 j = 1.589291662. 

This can then be used in equations (12. 2p and (12. 3p . 

Repeated use of (12. lip yields the equivalent series form of (12. 101) when q is an integer 



/3,E(-i) fc ( 

h — Cl \ 



2^ ^; - 1 



(3.4) 

fc=0 

When g is even this simplifies to 

(-) E(:)p)^f(-«M. 

k=0 V 7 V 7 fc=0 

When p = 3 and q = 3/2 the integrals in (I2.10p can be evaluated in terms of elementary 
functions but this leads to a transcendental equation for 03/2- Similarly when q is a half 
integer. 

As can be seen from (13.31) the numbers f3 q are not necessarily simple functions of q. For 
cases other than p = 1, 2, oo they can be numerically approximated and for this there are 
many other ways (I2.10p can be rewritten. For example, using a linear change of variables 
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and then the identities [16j 15.6.1, 15.9.21] the integral in (I2.10p can be written in terms 
of hypergeometric and associated Legendre functions. The result is 



2 _ U q-1 



dx = 2 q -\(3 q - l) q / [1 - (1 - f3 q )x/2} q ~ 1 x q - x dx 



[x 

2 q ~ l (B — l) q 
= ^ ^1(1 - g, g; 1 + g; (1 - /? fl )/2) 

= 2^r(g)(/3 g 2 -l) 9/2 p-\(/3 ? ). 

Numerical equation solvers can now be applied to any of these representations of 
jf 9 (x 2 - l)^ 1 dx in order to solve (l2~10l . 
The range of (3 q is known. 

Proposition 3.1. (5 q is a decreasing function of q and \/2 < /3 g < 2 for 1 < q < oo. 

Proof. A change of variables shows that ( 12.11) is equivalent to ff 9 1 y q ~ l {l + y)^ 1 ^ 2 dy = 
J(q), where J(q) = £ y*-\l - y)' 1 / 2 dy. Since (1 + y)' 1 / 2 < (1-y)- 1 / 2 for all < y < 1 
we must have (3 2 — 1 > 1. 

Equation (12.11) can also be written as I(q) + 2 f^(% 2 — l) 9 " 1 dx = J(q) where /(g) = 

Jq y' ?_1 (l + y) _1 / 2 dy. The argument above shows that J'(q) < I'(q) < for all 1 < q < oo. 
Since /3 q > \pT it must be a decreasing function of q. □ 

Even without knowing the exact value of (3 q we can use the method of Theorem 12.11 to 
obtain good estimates of the error in corrected trapezoidal rules. 

Theorem 3.2. Let c be the midpoint of [a, b}. Let f : [a, b] — > R such that f is absolutely 
continuous and J" G L p ([a,b]) for some 1 < p < oo. (a) If 1 < p < 2 let <p{x) = 
(x — c) 2 — {b — a) 2 /8. Equation (II. ip gives the quadrature formula 

(3-6) jf fix) dx = If (a) + /(&)] + lf'(a) - f'(b)] + E(f), 

where the error satisfies \E(f)\ < ||/"||i(& — a) 2 /16. (b) If 2 < p < oo let <fi{x) = (x — 
c) 2 — {b — a) 2 /12. Equation (II. ip gives the quadrature formula 

(3-7) f f{x) dx = if (a) + /(&)] + [f(a) - f(b)\ + E(f), 

where the error satisfies \E(f)\ < ||/"|| 2 (& - a) 2 ' 5 / '(12\/W) . (c) If 1 < p < oo let <p(x) = 
(x — c) 2 and 1/p + 1/q = 1. Equation (II. ip gives the quadrature formula 

(3.8) f f{x) dx = { -^- [/(a) + fib)} + [/'(a) - /'(6)] + E{f), 



2 L " v ' - WJ 8 
where the error satisfies 

\\f"Ub-a) 2 ^l q 



(3-9) < 



23+l/9(g+ 1/2)V9 ' 



s 
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Proof. Use the fact that L s ([a,b}) C L r ([a,b}) whenever 1 < r < s < oo. In (a) use the 
approximation from Corollary 12.21 and in (b) use the approximation from Corollary 12.31 
In (c) take a = 0, (3 = oo and then compute 



It is also possible to compute \\4>\\ q when (3 = 1 but this gives the trapezoidal rule (II. 2ft . 
Note that the coefficient in part (c) is strictly less than the corresponding coefficient in 
(ll.2j) for all 1 < q < oo. In the limit q — > oo the coefficient becomes 1/8 as in ( 11. 2ft . 



Let V m be the set of monic polynomials of degree m G N with real coefficients. Define 
F q :Vm — >■ M by F q (<f)) = \\<f)\\ q where 1 < q < oo and the norms are over compact interval 
[a, b]. Since F q (cf)) is bounded below for <fi G V m it has an infimum over V m . It also has a 
unique minimum at a polynomial that has m roots in [a, b]. As well, the error-minimizing 
polynomial is even or odd about the midpoint of [a, b] as m is even or odd. 

Lemma 4.1. (a) For m > 2, let <fi G V m with a non-real root. There exists if) G V m with 
a real root such that F q (if>) < F q (<f>). (b) Let <f> G V m with a root t ^ [a, b}. There exists 
ifi G V m with a root in [a,b] such that F q (ip) < F q (<f>). (c) If <f> minimizes F q then it has 
m simple zeros in [a,b\. (d) If F q has a minimum in V m it is unique, (e) F q attains its 
minimum over V m . (f) If (f) G V m is neither even nor odd about c := (a + b)/2 then there 
is a polynomial ip £ 'Pm that is either even or odd about c such that F q (ip) < F q ((f>). (g) 
The minimum of F q occurs at a polynomial <f> G V m with m simple zeros in [a, b}. If m is 
even about c then so is <f>. If m is odd about c then so is <f>. This minimizing polynomial 
is unique. 

Theorem 12.11 and its corollaries used only V% but the lemmas give results for V m for all 
m. These will be useful for considering integration schemes generated by polynomials of 
degree m, which we do not include here. 

The results of Lemma I4TT1 are well known and go back to Chebyshev and Fejer. To keep 
the paper self contained we have provided elementary proofs. For a full exposition and 
references to the original literature see, for example, j5], [7] and [15]. Three cases of the 
minimizing problem of Lemma 14.11 appear in the literature. When q = oo the minimizing 
polynomial in Vi is the Chebyshev polynomial of the first kind <f)(x) = x 2 — 1/2 = T 2 (x) /2. 
When q = 2 it is given as a Legendre polynomial <j)(x) = x 2 — 1/3 = (2/3)P2(x). When 
q = 1 it is the Chebyshev polynomial of the second kind <f>(x) = x 2 — 1/4 = U2(x)/A. These 
are for the interval [—1,1]. A linear change of variables is used for other intervals. These 
are all types of Gegenbauer polynomials, which are orthogonal on [—1, 1] with respect 
to a certain weight function. Gillis and Lewis [TT] give an argument to show that the 
minimizing polynomials are most likely not orthogonal polynomials for other values of q. 

Proof, (a) Write (f>(x) = [(x — r) 2 + s 2 ]u(x) for some r, s G E, s ^ 0, and uj G V m -2- Let 
if)(x) = (x — r) 2 u(x). For all x G M such that cj(x) ^ we have 1^(^)1 < 10(^)1- Hence, 




The rest follows as in the proof of Theorem 12.11 



□ 



4. Lemmas 



U\U<U\U- 
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(b) If m — 1 let 4>{x) = x — t. Direct calculation shows the unique minimum of 
q occurs at t — c. If m > 2, by (a) we can assume all the roots of <f> are real and 

write (f>(x) — (x — t)u(x) for some t £ [a,b], and u) G V m -\. Suppose t < a. Let 
if)(x) = (x — a)co(x). For all x G [a, b] such that co(x) ^ we have 1^(^)1 < 10(^)1- Hence, 
||V||,< U\\ q - Similarly if t> b. 

(c) Consider if){x) = (x — t) 2 with t G (a, b) and if) e (x) = (x — t + e)(x — t — e) = if)(x) — e 2 . 
For all x we have 1^(^)1 < 1^(^)1 + ^ an d for x ^ (t — e, t + e) we have 1^(^)1 < ^(^Ol- 
This shows that Halloo < Halloo if e > is small enough. Factoring now shows the zeros 
of any minimizing polynomial are simple. Similarly for t = a or b. 

For q = I, 

\\ipeh = / [e 2 - ip{x)] dx + / [if>(x) - e 2 ] dx 

< ||^|| 1+4e 3_ e 2 (6 _ a) 

< 111 f° r small enough e > 0. 

And, if g G L°°([a, b}) such that \g\ > almost everywhere then ||^ e <3>||i < H^Hi for small 
enough e > 0. A similar construction is used when t equals a or b. Factoring if) shows the 
zeros of any minimizing polynomial must be simple. 

For 1 < q < oo use the same construction, with Taylor's theorem in ( 14.11) . to get 

Ue\\ q g = [ [e 2 - iP(x)] g dx + [ [i){x) - e 2 ] q dx 

2? An. _L / ,IA(rA - \nl,(rA - . 2 1 ^ 



(4.1) < / e 2q dx + if) q (x) - qe 2 [if)(x) - e 2 ] 9 dx 

J xe(t—e,t+e) Jxg(t—e,t+e) 

< 2e 2q+1 + \\iP\\ q q - qe 2 f (Se 2 )^ 1 dx 

J x<^{t-2e,t+2e) 

= 2e 2q+1 + \\^\\ q [-3 q - 1 qe 2q {b-a-4e) 

< \\ip\\g for small enough e > 0. 

As with the q = 1 case above, it now follows that zeros of minimizing polynomials must 
be simple. 

(d) Suppose the minimum of F q occurs at both 0, u G V m . Let if) = ((f) + u))/2. Then 
if) G V m so \\(f>\\ q < \\if)\\ q < (\\<P\\g + IMIg)/2 = ||0|| g - The Minkowski inequality must 
then reduce to equality. For 1 < q < oo this means <f> = du for some d > [131 p. 47]. 
Since 0, u) G P m we must have d = 1. If g = 1 then there is equality in the Minkowski 
inequality if and only if (f>u) > 0. If there is equality, then <f> and if) must share roots of odd 
multiplicity. But by (c), <f> and if) have m simple zeros in [a, b}. Hence, <f> = if). For q = oo, 
Lemma I4.2I shows there are a < x\ < ■ ■ ■ < x m _i < b, a < y\ < ■ ■ ■ < y m -\ < b and 
a < z\ < ■■ ■ < z m _ x < b such that \<p{xi)\ = \u)(yi)\ = \ip(zi)\ = \\<f)\\oo, (f>'( x i) = = 
if)'( Zi ) = for each 1 < i < m - 1. Let = {x^ 1 , M u = {y l }T = ~ 1 \ = {z^ 1 . 
Let z G then ll^ll^ = \if)(z)\ = \<f>{z) + w(z)\/2 < ||0||oo unless z G (1 M u . Hence, 

= M^- Therefore, <f>'(xi) = if)'(xi) for each 1 < % < m — 1. And then <f> = u). 

(e) By parts (a) and (b) we need only consider <f> G V m with m real roots in [a, b]. Let 
U G [a, 6] for 1 < % < m and define (f> G P m by <f>(x) = YYh=i( x ~ U)- 
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Suppose 1 < q < oo. Define G q (t) = J^YliLi\ x ~ U\ q dx. Then G q : [a,b] m — > R attains 
its minimum over [a, b] m if and only if F q attains its minimum over V m - The set [a, b] m 
is compact in R m . And, G q is continuous. For, suppose t G [a, b] m and s^ fc ^ G [a, 6] m for 
each k G N such that s^ fc ^ — >• t in the Euclidean norm. We have 

\\\x- sf ] \ q dx< / ]J(b-a) q dx= (b-a) mq+1 . 
- i=i ^ a i=i 

By dominated convergence (for example, [TJ Theorem 7.2]), lim^oo G q (s^) = G q (t) and 
G 9 is continuous. Therefore, G q attains its minimum over [a,6] m . 
The case q = oo is similar, using F^cj)) = max x . e [ a 6 ] niiil x ~U\- 

(f) Without loss of generality, b = —a. Suppose ip G V m is the unique minimizer of 
F q . Let u(x) = ip(—x) if m is even and u(x) = —ip(—x) if m is odd. Then u G P m and 
ll^ll? = llV'llg- Let ( = (ip + u)/2. Then £ G ? m and is even if m is even, odd if m is odd. 
Also, \\C\\q < (HV'llg + IMIg)/2 = HV'llg- Hence, ( = ip = u. But then ip is even or odd as 
m is even or odd. If G V m is neither even nor odd then F q (ip) < F q ((p). □ 

Lemma 4.2. Suppose <p G V rn is a minimum of F^. Then <p{x) = YliLi( x ~ U) f or 
a < ti < t 2 < ■ ■ ■ < t m < b. For each 1 < i < m — 1 there is ^ G (tj, s^c/i i/iai 

= ii^iioo. 

Proof. This follows from the q = 00 case of Lemma [4.1( c). □ 



5. /" Henstock-Kurzweil integrable 

Let iJfr([a, 6]) be the set of functions integrable in the Henstock-Kurzweil sense on [a, b}. 
See, for example, [T2J. Note that L s ([a, b}) C L r ([a, 6]) C M([a, 6]) for all 1 < r < s < c». 
An example of a function / G HK"([0, 1]) that is not in any L P ([0,1]) space is / = F' 
where F(x) = x 2 sin(a;~ 3 ) for x G (0, 1] and F(x) = 0. In this section we use the method 
of Theorem 12.11 to choose a monic quadratic <ft to minimize the error in the resulting 
corrected trapezoidal rule when /" G HK([a, b]). We also consider the case when /' is 
merely continuous and then /" exists as a distribution. For all of these cases, it turns out 
that amongst corrected trapezoidal rules (jl.ip . the trapezoidal rule itself minimizes the 
error. 

If / G HK([a,b]) then the Alexiewicz norm of / is ||/|| = sup x6 r a6 ]|J a f(t) dt\. If 
g : [a, b] — > R is of bounded variation then fg G HK([a,b]). Integration by parts shows 
that I J f(x)g(x)dx\ < \J f\\g(b)\ + ||/||Vgr, where Vg is the variation of g. This is a 
version of an inequality known in the literature as the Darst-Pollard-Beesack inequality. 
See [6], [2]. However, it appears earlier in |18j . It is proved for a more symmetric version 
of the Alexiewicz norm for Henstock-Kurzweil integrals in [191 Lemma 24] . 

Under the Alexiewicz norm, HK([a,b}) is a normed linear space but is not complete. 
We will discuss integration in the completion later in this section. 

Theorem 5.1. Let c be the midpoint of [a,b\. Let f :[a,b] — > R such that f G C([a,b}) 
and f" G HK([a, b]). Amongst all monic quadratic polynomials used to generate (11.11) . 
taking <j>(x) = (x — c) 2 — (6 — a) 2 /4 = (x — a)(x — b) minimizes the error \E(f)\. This 
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gives the quadrature formula 

a 



(5-1) J f(x)dx = —[f(a) + f(b)] + E(f), 

where 

(5-2) \E(f)\< 11/1(6 -af/A. 

For a uniform partition given by X{ = a + (b — a)i/n, < % < n, the composite formula is 
the trapezoidal rule 



r>b 

(5.3) / f(x)dx 



n-1 



wrw-a) 2 



+ E(f), 



2n 

where 

w)i * 4n ■ 

The coefficient of \\f"\\ in the error bounds is the best possible. 

Proof. From flTT) we have \E(f)\ < (\f*f"(x) dx\\(j)(b)\ + \\f\\V<f>)/2. Write <f>(x) = 

(x - r) 2 + s. Then V(f) = 2 £\ x — r\dx. This is minimized when r = c. Note that 
4>(b) = if and only if s = — (b — a) 2 /4. The unique error-minimizing polynomial is 
then <j)(x) = (x — c) 2 — (b — a) 2 /4 — (x — a)(x — b), for which V<p = (b — a) 2 /2. Hence, 
|^(/)| — ll/"ll(^ — a) 2 /A. To show this estimate cannot be improved, let (ip n ) be a delta 
sequence as in the proof of Corollary [2721 Define f"(x) = i J n(x — a) — 2ip n (x — c)+ip n (b—x). 
Then f' n \x) dx = and Wf^W = 1- Integrate to get f n (x) = J^(x — t)f£(t) dt, modulo a 
linear function. Since <\> is continuous, lim^oo | f"(x)4>(x) dx\ = \4>{a) — 20(c) + 4>(b)\ = 
(b-a) 2 /2. "' □ 

In the terminology of Theorem (12.11) the minimizing polynomial has (3 = 1 and a = 
(b — a)/2 and is the same as the one that generated the trapezoidal rule in (jl.2p . Notice 
that the variation is additive over disjoint intervals so the error in the composite formula 
is order 1/n, compared with error of order 1/n 2 in Corollary 12.51 when /" 6 L p ([a,b}). 

An equivalent norm is ||/||* = sup/| f T f(x) dx\ where the supremum is taken over all 
intervals I C [a, 6]. It is easy to see that ||/|| < ||/||* < 2||/|| for all / E HK([a,b\). In 
terms of this norm, (15. 2 j) implies \E(f)\ < — a) 2 /A. This inequality appears as 

Theorem 21 in (TO]- We improve this inequality by a factor of 1/2 and show our result is 
sharp. 

Corollary 5.2. With the assumptions of Theorem \5.1\ we have \E(f)\ < (6 — a) 2 ||/"||*/8. 
The constant 1/8 is the best possible. 

Proof. Use the second mean value theorem for integrals [3j. If is monotonic on [a, b] 
then there is £ G [a, b] such that 

E{f) = \ f f"(x)<Kx)dx = ^- f f'\x)dx + ^ f f"{x)dx. 
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Then \E(f)\ < (\<f>(a)\ + |0(&)|)||/"||72. This is minimized by taking <f>(x) = (x - a) 2 or 
<f)(x) = (x — b) 2 , for which \E(f)\ < (b — a) 2 ||/"||*/2. If <ft has a minimum at r G (a, b) 
then there are £1 G [a, r] and £ 2 G [r, b] such that 

E (f) = Mj^ r{t) ^ + m£ nx) ^ + mj^ r(t) ^ + m^ r{x)dx . 

It follows that < (|0(a)| + |0(r)| + |0(6)|)||/"||*/2. Choosing <f>(x) = (x - a) (a - 6) 

minimizes the coefficient of ||/"||* and we get \E(f)\ < {b — a) 2 ||/"||*/8. To prove this 
estimate is sharp, let (ip n ) be a delta sequence as in the proof of Corollary 12.21 Define 
f;{x) = V„(x - c). Then \\f^\\* = 1 and \E{ip n )\ -> (c - a) (6 - c)/2 = (6 - a) 2 /8. □ 

The completion of HK([a,b}) in the Alexiewicz norm is the Banach space *4. c ([a,o]). 
Each element of *4 c ([a, b]) is the distributional derivative of a function in B c ([a,b\) = 
{E G C([a, b]) | F(a) = 0}. Note that B c ([a, b]) is a Banach space under usual pointwise 
operations and the uniform norm. If / G «4 c ([a, 6]) then there is a unique primitive 
F G B c ([a, b]) such that F' = f. The distributional derivative is (F',4>) = —{F,<f>') = 
— J a F(x)4>'(x) dx where <p £ C£°(l£)- The Alexiewicz norm of / is ||/|| = sup^gu wl/^ f\ = 
H-Flloo- This makes ^4 c ([a, 6]) into a Banach space isometrically isomorphic to B c ([a,b]). 
The continuous primitive integral of / G *A c ([a, &]) is then f^f = F(b) — F(a). If g 
is of bounded variation then the integration by parts formula is given in terms of a 
Riemann-Stieltjes integral as J^fg = F(b)g(b) — J F(x) dg(x). Note that ^4 c ([a, 6]) 
contains HK([a, b}) and hence L p ([a, b}) for each 1 < p < oo. If F is a continuous function 
such that its pointwise derivative F'(x) = almost everywhere then the Lebesgue integral 
of F'{x) exists and is but the continuous primitive integral is j a F' = F(b) — F(a). If F 
is continuous such that the pointwise derivative exists nowhere then the Lebesgue integral 
of F'{x) does not exist but F' G A c ([a,b}) and f*F' = F(b) - F(a). See [20] for details. 

Note that if F G C([a,b}) then F' G A c ([a,b]) and f* F' = F(b) - F(a). An advantage 
of the continuous primitive integral is that the space of primitives is simple. It is the 
continuous functions while for the Henstock-Kurzweil integral it is a complicated space 
called ACG*. See [5] for the definition. 

Corollary 5.3. Let f G C 1 ([a, b]). Then f" G *4 c ([a, b)) and the formulas in Theorem \5.1\ 
hold. 

In a sense this now reduces to estimates on /' since the Alexiewicz norm of /" is the 
uniform norm of /'. The formulas in Theorem 15 . 1 1 also hold when /' is a regulated function. 
This is a function that has a left limit and a right limit at each point. See [21] for details. 

6. Exact for cubics 

In this section we show ( 11. II) is exact for all when / is a linear function. We also show 
(11.11) is exact for all cubic polynomials / whenever <p{x) — (x — c) 2 — (6 — a) 2 /12. 

Theorem 6.1. Let f : [a, b] — > K and let c be the midpoint of [a,b]. Let (ft be a monic 
quadratic. Write 

(6.1) f b f(x) dx= l - [-/(a)0'(a) + f{b)<f/{b) + /'(a)0(a) - f{b)<j>{b)\ + E{f). 
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(a) If f is a linear function then E(f) = for all such 0. (b) E(f) = for all cubic 
polynomials f if and only if <j)(x) = (x — c) 2 — (b — a) 2 /12. 

Proof. The proof of (a) is straightforward using E(f) = (1/2) J b f"(x)<j)(x) dx. By (a) 
we need only consider f(x) = Ax 3 + Bx 2 . Write <p(x) = x 2 + Cx + D. The equation 
J^(6Ax + 2B)(x 2 + Cx + D) dx = gives the linear system 

(6.2) 2{a 2 + ab + b 2 )C + 3{a + b)D = --(a + b){a 2 + b 2 ) 

(6.3) (a + b)C + 2D = -\(a 2 + ab + b 2 ). 

The solution is C = — (a + b), D = (a 2 + Aab + b 2 )/6, and this gives <f)(x) = (x — c) 2 — (b — 
a) 2 /12. ' 1 □ 

Note that this is the same <fi as in Corollary 12.31 Also, E(f) = for all quadratic 
polynomials / if and only if (16.31) holds. 
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